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Abstract. We study linear power spectra and formation 
of large scale structures in flat cosmological models with 
A > and cold plus hot dark matter components. We re- 
fer to these models as mixed A models (MLM). The hot 
component consists of massive neutrinos with cosmologi- 
cal density Qh and number of neutrino species as a free pa- 
rameter. The linearized Einstein-Boltzmann equations for 
the evolution of the metric and density perturbations are 
integrated for a set of values of the cosmological parame- 
ters. We study MLM models with present matter density 
in the range 0.25 < Om < 1, dimensionless Hubble con- 
stant 0.4 < h < 0.7 and the hot dark matter content with 
a ratio within the limits < Hh/^m < 0.3. For all the 
considered models we assume a scale-invariant primeval 
spectrum. 

The density weighted final linear power spectra are 
normalized to the four year COBE data and have been 
used to constrain the parameter space by a comparison 
of linear predictions with the current observational data 
on large scales. The consistency of MLM predictions with 
the observable data set is best obtained for models with 
one species of massive neutrinos and Qh/Qm < 0.2. Of 
the considered linear tests the strongest constraints on 
Qm that we obtain arise by comparing the cluster X-ray 
temperature function with that observed at the present 
epoch. Consistency with the estimated cluster abundance 
can be achieved for COBE normalized MLM models with 
Qh/Om < 0.2 and 0.45 < Om < 0.75 for h = 0.5. If 
h = 0.7 then 0.3 < Qm < 0.5. These constraints are at ler 
level and standard MDM models are clearly ruled out. 

We note that the range of allowed values for SIm, that 
we obtain for MLM models from linear analysis, is also 
approximately the same range that is needed in order to 
consistently satisfy a variety of independent observational 
constraints. 
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1. Introduction 

In the standard framework of gravitational instability the- 
ory present day structures must have been formed through 
the growth of small inhomogeneities from an initial ran- 
dom Gaussian density field present at very early epochs, 
with a scale- invariant Harrison-Zel'dovich spectrum. Clus- 
tering analysis of the large scale structure in the Universe 
has been improved in recent years by observations of the 
spatial distribution of galaxies and cluster of galaxies, as 
well as cosmic microwave background (CMB) anisotropies. 
Thus any theory that wants to fit the observed large scale 
clustering must be consistent with a set of constraints 
over more than three decades in length: from galaxy cor- 
relation (~ lh~ 1 Mpc, H = 100hKmsec^ 1 Mpc~ 1 ) up 
to the quadrupole CMB anisotropies detected by COBE 
(~ SOOOh^Mpc, |5moot ct alT" 



1992) 



In the standard FRW metric the fundamental back- 
ground cosmological parameters are related by a mutual 
relation. It is also understood that for these parame- 
ters (the present matter density Qmi the value of the 
Hubble constant Ho and the age of the universe to) the 
range of values allowed by observations must be consis- 
tent with their FRW relation. If we adopt the standard 
inflation theory we shall assume that the total energy 
density in the Universe is equal to the critical density 
(Oo = !)■ On the other hand, observations and theoreti- 
cal predictions suggest that the amou nt of baryon density 
must be small (Q b ~ 0.05(fe/0 .5)~ 2 , |Walker et al. 1991 
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Copi, Schramm fc Turner 1995 ). A much larger contribu- 
tion to the matter density must be of non-baryonic nature. 
This is so-called dark matter (DM) problem and it has 
observational support from dynamical estimates (SIm ~ 
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0.2 — 0.3, Bahcall 199E and references cited therein). It 
must be stressed that there is not yet a firm evidence for 
f2jvf = 1, this value for Qm being required by the inflation- 
ary paradigm. The simplest possible model of DM is the 
one in which the universe is dynamically dominated by 
a single massive collisionless particle. The most popular 
DM candidates are collisionless massive particles, which 
decoupled from cosmological plasma either when being rel- 
ativistic (DM particles like massive neutrinos - Hot DM) 
or non-relativistic (hypothetical massive particles - Cold 
DM). 

The simplest model is CDM, where the power spec- 
trum of the perturbations depends on a single param- 
eter, the present cold density flc- Historically the first 
model to be considered was a neutrino with a non-zero 



scales ( Klypin ct al. 1993f). For the c luster correlation 



rest mass HDM ( [Bond fe Szalay 1983^ |Zakharov 1979D 
The HDM model was soon rejected because nu 
merical 



Davis 1983) pro- 
in agree- 
model 



simulations ( White, Frenk 
duced nonlinear structures too late to be 
ment with QSO existence. The standard CDM 
(SCDM) has been analysed in considerable details 
( Davis ct al. 1985] ; Davis 
cited therein). It 



Efstathiou 1988 



and refer- 
galaxy cluster- 



ences cited therein), it can match 
ing (^5 10h~ 1 Mpc) with a bias parameter b g lZ 2, al- 
though it lacks of sufficient power on large scales. 
The SCDM model however became seriously inconsis- 
tent with clustering data after the COBE detection of a 



quadrupole anisotropy in the CMB (Smoot et al. 1992) 
When the power spectrum of the model is normal 
ized to the COBE data SCDM models have un 



avoidable difficult i es (White, Efstathiou fe Frenk 1993a ; 
Olivier et al. 1993[ |Mo, Jing fe Borner 1993j ; Jing fe Val- 
darnini 1993; Dalton et al. 1991 ; Baugh & Efstathiou 
1994). 

The Mixed Dark Matter models (MDM) was at 
first proposed and discussed in a few papers (Fang, Xi- 
ang & Li 1984; Shafi & Stecker 1984; Achilli, Oc- 
chionero & Scaramella 1985; Valdarnini & Bonometto 
1985) as an example to overcome the standard Hot Dark 
Matter (SHDM). Late r SMDM models were readr essed 
again ( [Holtzman 1989| ; |Van Dalen fc Schaefer 1992| ). Fi- 
nally, it was proposed that a certain mixture pro- 
portion of about 30/70% for the Hot/Cold DM 
QSchaefer fc Shafi 1992|;pavis, Summers fc Schlcgcl 1992] ; 



Taylor fc Rowan-Robinson 1992] ; Klypin ct aT 1993; 
Pogosyan & Starobinsky 1993) as a DM model satisfy- 
ing a wide body of observational data for clustering ( dis- 
tributions of galaxies, galaxy clusters, quasars, Ly-a sys- 
tems, etc.) on typical scales < Ils ~ 100 — 150h~ 1 Mpc 
and with an Harrison- Zeldovich spectrum of the primor- 
dial cosmological perturbations on scales > Ils- Be- 
cause of the changes introduced into the power spec- 
trum by neutrino clustering, for SMDM models it was 
shown the possibility to reconcile the evidence of high 
coherent velocity flows on large scales (~ 50h~ 1 Mpc ) 
with the moderate galactic pair velocities on megaparsec 



function ( Holtzman & Primack 1993: Jing ct al. 1993 



Klypin fc Rhce 1993[ ) the two-point function for SMDM 
models appears to be consistent with data for R^, 20 — 
30/i _1 Afpc ( on the contrary to what found for CDM ). 

Analytical approximations for the present epoch trans- 
fer functions are given by Holtzmann (1989); Pogosyan & 
Starobinsky (1995); Ma (1996); Eisenstein & Hu (1997). 
Early numerical simulations have considered £Ih = 0.3, 
we will use the notation H for massive neutrinos, but 
this model does not produce Ly-a systems as much as 



observed (Mo & Miralda-Escude 1994). For this reason a 
value of CIh^ 0.2 results in a better fit (Klypin et al. 1995; 
Ma & Bertschinger 1994). In comparison with SCDM or 
SHDM, MDM models are more complicated, we have the 
choice of two independent parameters : the ratio Qh/^c 
and ma ( if one allows for more than one species of massive 
neutrinos). The spectral index n of the post-inflationary 
spectrum is taken n = 1. It is possible to consider also the 
role of gravitational waves, which can change significally 
the normalization of the spectrum and the formation of 
structures (e.g., Ma 1996 ). 

It must be stressed that SMDM models are in dif- 
ficulty with the present upper limits on the age of 
the Universe: if we assume h = 0.5 and the age of 
globular clusters greater than IbGyr, in the case of a 
flat FRW model the needed value of is < 0.6. 
The other difficulties are connected with later galaxy 
and quasar format ion ( Pogosyan fc Starobinsky 1995 ; 



Cen fc Ostriker 1994! ) . Another difficulty for SMDM mod- 
els is that when the power spectra are normalized to the 
COBE 4-yr data the linear theory overpredicits the ob- 
served cluster abundances. The estimated uncertainties in 
the normalization and linear calculations can hardly fit 
the cluster number density within present data error bars. 
The difficulty can be reduced if one removes the constraint 
of a scale invariant spectra and introduces a small tilt ( 
n~ 0.8 -0.9, |Ma 1996j ). 

In alternative to MDM models spatially flat low den- 
sity models with a positive cosmological constant are 
considered to be a viable generalization of SCDM after 
COBE. These models are termed ACDM and were consid- 



ered even before the quadrupole detection ( |Pccblcs 1984 
and references therein). There are two main reason to 
consider ACDM models as an alternative to MDM: the 
present lower limit on the age of the Universe and the 
baryon fraction in galaxy clusters. 

Current uncertainties for the age of globular clusters 
give to = 15 ± 2Gyr (Chaboyer et al. 1996 ). If one takes 
a lower limits to > ISGyr then in the fijvf = 1 case 
the maximum value for h is 0.5. This value is below the 
ran ge allowed by recent HST measurement : h — 0.7 ± 
.1 ( Freedman et al. 1994 ; Reiss, Kirshner fc Press 1995 ). 
The introduction of a positive cosmological constant al- 
lows h to be higher for a fixed age to than in CIm = L 
The other reason to consider ACDM model is from X- 
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ray observations of galaxy clusters. If clusters are a rep- 
resentative sample of the matter content in the Uni- 
verse, then their baryons/total matter ratio should not be 
different from what expected from standard nucleosyn- 
thesis QWhitc et al. 1993b| ). For h> 0.5 this is achieved 
if Hm~ 0.35 ± 0.2. Quite interestingly this range is 
close to that obtained from dynamical estimates. After 
COBE detection of the quadrupole anisotropy ACDM 
models have been further considered (Kofman, Gnedin 
& Bahcall 1993; Klypin, Primack & Holtzman 1996; 
Liddle et al. 1996a; Liddle et al. 1996b). They can 
fit several constrai nts on large scales (k^S lh~ 1 Mpc~ 1 , 
Liddle et al. 1996b), but small scales clustering it is excess 
(Klypin, Primack fc Holtzman 1996) by a factor ss 2 — 3 



when compared to estimates from galaxy catalogs, so this 
is a difficulty of the model which can be solved with the 
introduction of an antibias of small scales. A possibility 
which seems unplausible ( Primack fc Klypin 1996| ). 

These difficulties have suggested that consistency with 
present data can be achieved for standard DM models with 
the introduction of one extra parameter. Possible variants 
are two species of massive neutrinos ( Primack 1997 ), or 
a tilt of the primordial spectrum ( Cole et al. 1997 ). 

The main aim of this paper is to consider an alternative 
possibility , that is MDM models with n = 1 and a non- 
zero cosmological constant (MLM) . These models have in 
fact the following free parameters : VLh/^Im, and the 
number of species of massive neutrinos. If one considers 
one species of massive neutrinos, then MLM are on the 
same foot of the previous alternatives, the main advantage 
being that they retain the inflation paradigm and a scale- 
invariant initial spectrum. 

In order to study these models we have numeri- 
cally integrated forward in time the linearized Einstein- 
Boltzmann equations for the evolution of the metric 
and density perturbations. The final power spectra have 
been normalized according to the four year COBE data 



(Bunn & White 1997). We have then used linear pertur- 



bation theory to find the parameters for which MLM 
models are consistent with a variety of current obser- 
vational constraints on large scales (power spectrum, 
cluster-cluster correlations, bulk velocities, cluster abun- 
dance). In order to make comparison with observations 
we will consider the following set of parameter space: 
0\ = 0,0.31,0.45,0.65,0.74; VL H /VL M = 0.1,0.2,0.3; num- 
ber of massive and massless neutrinos (3h — 2,4,6 and 
f3 v — 0,2,4,6 ( here (3 is the sum of spin states ). In 
some cases we considered a number of neutrinos species 
greater than three, this possibility being of no physi- 
cal relevance, we ran these models for studying how dif- 
ferent numbers of massive and massless particles would 
change the spectra evolution. Our final analysis will be 
restricted to one massive neutrino. The fifth parameter is 
h(h = 0.4,0.5,0.63,0.7). 

The outline of the paper is as follows: in Section 2 we 
will present our model description and equations describ- 



ing the linear perturbations in MLM models. In Section 3 
we will study the power spectra dependence on model pa- 
rameters using COBE 4-year data normalization. In Sec- 
tion 4 some linear tests are given for MLM models, in- 
cluding our results for the cluster mass and temperature 
function. Our conclusions are given in Section 5. 

2. Theoretical framework 

2.1. Equations and method 

We have studied the evolution of density perturbations 
in a flat Friedmann cosmology. In our models the present 
total density in critical units is ilo = &c + + = 
1 = + ^a- The notation H means massive neutrinos, 
C cold dark matter and Ha = A/(3Hq). 

In our integrations we have followed numerically the 
time evolution of the linearized Einstein-Boltzmann equa- 
tions for the metric and adiabatic density perturbations. 
We have considered scalar modes only and treated the per- 
turbations for the following particle species : cold, massive 
and massless neutrinos {v), photons. We treat baryons and 
radiation as a single ideal fluid. We have assumed that all 
collisionless particles were decoupled from radiation be- 
fore the beginning of our computation. In a collisionless 
medium the pressure anisotropy is different from zero and 
the equations for density contrast and velocity are not 
sufficient to describe the perturbations in H and v par- 
ticles. Collisionless components must be described by the 
Boltzmann-Vlasov equation. The density contrast, flux ve- 
locity and pressure anisotropy are given by the moments 
of the perturbed distribution functions. 

The set of equations for fluid and collisionless media, 
with the Einstein equations for metric perturbations, de- 
scribes the evolution of density perturbations in MLM 
models. The method of numerical solution of this system 
is described in Valdarnini & Bonometto (1985), and we 
refer to this paper for more details. The generalization of 
the equations to include a non-zero cosmological constant 
is straightforward. 

Our numerical calculation is done for the following 
range of perturbation masses : from M = 1O 2O M0 down to 
M = 10 11 M Q . The integrations start at the initial redshift 
Zi = 10 9 and stop at the final epoch Zf — 5. Because of 
the large region of parameter space that we have spanned 
with our integrations, and the computational resources 
we had available, we have decided to stop the numerical 
integrations at a final epoch Zf > 0. As a compromise be- 
tween computing budget and the accuracy needed to eval- 
uate observational linear variables at the present epoch we 
chose Zf = 5. 

The initial conditions for our linear computations are 
given by the Harrison- Zeldovich ( n — 1 ) spectrum of 
density fluctuations: 



Sc(k) 



A.i k . 



(i) 
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where k is the comoving wavenumber of the perturbation 
and At is an arbitrary constant. The total matter density 
perturbation is then defined as : 

6m = &PC t = J-(«c<fc + n H S H ). (2) 
Pc + Ph 



The transfer function T(k) can be defined as the ra- 
tio of the amplitude of the Fourier mode Sm(zj, k) to the 
one of minimal k, which corresponds to the wavenumber 
of the maximal mass perturbation M = 10 20 M Q . The fi- 
nal linear transfer functions of our integrations have then 
been evaluated at z = using analytical formula. The ap- 
proximation involved neglects the changes in the shape of 
the transfer functions that take place between z = 5 and 
z = because of the decrease in the neutrino streaming. 
The error involved in the computation of the linear vari- 
ables, with which we test the models against a set of data, 
is negligible in most of the considered cases and amounts 
to a few percent when Qh/Qm — 0.2 and Qm = 0.25. 

To describe the evolution of fluctuations between Zf 
and z = we have applied the well known exact analyti- 
cal solutions of Einstein equations for a perturbed flat dust 
model with non-zero A ( |Kofman fc Starobinsky 1985j 
Lahav et al. 199l[). 



A non-zero cosmological constant changes the expan- 
sion rate so that the scale factor has the following expres- 
sion 



a(t) 



n A 



1/3 



smh 2/3 (-V / Olffoi), (3) 



and as a result the growth of perturbations relative to that 
of a critical density universe, after the moment of equality 
of the cosmological constant density to the matter one, it 
is suppressed. The suppression coefficients for density and 
velocity perturbations are the following: 



K s (t) 



1 



a(t)dt 



(4) 



i * 5 / a a 

K vV) = o - — 



3 W ' ' ' <t)dL 
The total suppression at z=0 is well approximated by 



(5) 



1 

70 



209fi M - 



140 



Om 4/7 



, (6) 



for density perturbations flCarroll, Press fc Turner 1992 ) 



and by f^M ' 6 for peculiar velocities (Lahav et al. 1991) 



We express the final matter power spectrum which 
gives rise to the observed large scale structure of the Uni- 
verse as 

P(k) = AkT 2 (k)K 2 s (t )/n 2 Ml (7) 

where A is a normalization constant and T(k — ► 0) — > 
1 is the total matter transfer function. In Eq.(7) P{k) 
is the density weighted power spectrum, i.e. P(k) — 



ZlcPH 2 



2.2. Normalization procedure 

In order to constrain a cosmological model one has to 
fix the normalization amplitude of the density fluctuation 
power spectrum at the present epoch. 

The most accurate method for normalizing P(k) is to 
make use of the COBE satellite data for cosmic microwave 
background (CMB) anisotropies (3moot et al. 1992, 



Bennett et al. 1994[ |Bennctt et al. 199J ). Bunn & White 
(1997) have carried out a likelihood analysis of the 4-year 
COBE DMR sky maps and have obtained the best-fitting 
quadrupole, the value of which depends strongly on the 
initial index n, Hm and the space curvature. For example, 
for a Sachs- Wolfe n = 1 spectrum < Q >= 18.7 ± 1.3(J,K, 
for n — 1.5 < Q >— 13 ± 0.8/j.K; in the case of flat cosmo- 
logical models with non-zero A and open models its value 
increases when VLm decreases. For the normalization of 
the spectra in different cosmologies Liddle & Lyth (1993) 
have proposed instead of the quadrupole to use the ampli- 
tude of the density perturbation at horizon-crossing 5hor 
defined by 



A 2 (fc) 



k 3 P(k) 



2tt 2 



= 5 



hor 



ck 



T 2 (fc), (8) 



which fixes the present-day normalization of spectra. A fit 
to the 4-year COBE data for flat models with the n = 
1 post-inflation spectrum analysed here has the follow- 



ing simple form ( 


Liddle et al. 1996a 




Liddle & Lyth 1993 


Liddle et al. 1996b[ Bunn & White 


1997 


): 



<W("m) = 1-94 io-^ L 



5n-0.785-0.05[nSl u 



(9) 



Proceeding from the definitions of Shor and the power 
spectrum, Eq.(7), the normalization constant is calculated 
as 

c \ ( 
H ) \K S 



27r2 $hor 



(10) 



3. Power spectra 



The power spectra of Eq.(7), normalized to the COBE 4- 
year data, are depicted in Figs. 1-3 for different models. 

Fig. 1 shows how the number of species of massive 
neutrinos, when is kept fixed (= 0.2), modifies the 
amplitude of the spectra at k > O.lhMpc -1 . As we can 
see, increasing the number of species of massive neutrinos 
from 1 to 3 decreases the power at scale < 100ft. -1 Mpc 
of about 1.6 times. Increasing the number of species of 
massless neutrinos when the number of the massive ones 
is fixed also suppresses the power. Both of these effects 
have a simply explanation. 

The dynamics of the perturbations in the hot compo- 
nents depends on two characteristic scales: one of them 
is the free-streaming scale, which is about the inverse of 
the Hubble scale, when the particles first become non- 
relativistic. In our case fc# = kn(zH) — H(zh)/c ~ 
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K [ h Mpc-' ] K [ h Mpc-i ] 



Fig. 1. The linear density power spectrum P(k), normalized 
to the 4-yr COBE data, is shown for models with Qa = 0, 
Q.h/Qm = 0.2, h — 0.5 and different values of the (3h and /3„ 
parameters. The values chosen are: Ph = 2 f3 v = 4 (solid line), 
Ph = 6 fi v = (dotted line), (3 H = 4 /?„ = 2 (short-dashed 
line), /3ff = 2 /3„ = 6 (long dashed line) 



Fig. 2. We show the power spectra dependence on the Qh/^m 
parameter. The models chosen are MDM models with = 0, 
h — 0.5, (3h = 2, p v — 4 and different values of Q,h : 
fin =0.1 (solid line), Q H = 0.2 (dotted line) and fi H = 0.3 
(short-dashed line) 



Mpc^ 1 where zr is the redshift at which the H- 
particles become non-relativistic, run is in eV and H{z) 
is the Hubble function. 

The second one is connected to the free-streaming 
scale, i.e. the minimal size ~ fc#(z) -1 of if-particle ob- 
jects that can collapse at a given redshift z. After the 
epoch zh the velocities of massive neutrinos are redshifted 
away adiabatically, then kJ I l {z) decreases. On large scales, 
k < kn, gravitational instability develops in a standard 
way. For the growing mode of adiabatic perturbations 
S ~ Sh — Sc ~ a- On the other hand the growth of 
density perturbations on smaller scales, k > kn, is re- 
duced from the moment when a perturbation of a given 
scale enters the horizon and up to the time when the free- 
streaming scale becomes smaller than the perturbation 
size. Because of the free-streaming effects density pertur- 
bations are erased away in the part of <~ Hh, they still 
grow in the cold component in such a way that the total 
density growth, although less than in CDM , appears much 
larger than in HDM model. Finally for k < kn{z) the 
perturbations in both components develop similarly: the 
dynamical amplitude 8r gradually approaches the other 
component Sc- 

Thus decreasing the mass of massive neutrinos (in- 
creasing their number of species when is fixed) shifts to 
later epochs the time when they become non-relativistic 



T 1 1 1 — I ' ' ' J_ 1 1 1 1 — I I I I | 1 1 1 1 — I I I I- 




K [ h Mpc" 1 ] 



Fig. 3. Power spectra are shown for MLM models with h = 0.5, 
/3h = 2,/3„ = 4, Q.h/Qm = 0.2 and different values of the 
f?A parameter: Qa = 0. (solid line), S7a = 0.31 (dotted line), 
Qa = 0.45 (short dashed line), Qa = 0.65 (long dashed line), 
Qa = 0.74 (dot-short dashed line) 
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and the collisionlcss damping becomes more effective in 
slowing down the growth of matter density perturbations. 

Increasing the number of species (N v ) of massless 
neutrinos results in a longer duration of the radiation 
dominated era and in a smaller redshift a~ q = 1. + 
z eq ~ 4.16 l0 4 n M h 2 /(1+0.227JV,,) of the matter-radiation 
equality. It is useful to define the comoving wavenumber 
k with respect to the scale that crosses the horizon at the 
matter-radiation epoch, i.e. 

q = k/k eq = k/[QAn M h 2 {\ + 0.2277V I/ )- 1 Mpc _1 ]. 

Because of the approximate scaling of the transfer function 
with the dimcnsionlcss parameter q, an increase in N v 
implies a reduced amplitude of the power spectrum at a 
given k. 

So both of these effects, increasing the number of mas- 
sive and massless neutrinos, increase the duration of the 
epoch when the growth of density perturbations is less 
effective. The dependence of the spectra on Qh, or on 
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Fig. 4. Here P(k) is plotted for two groups of MLM models 
with fixed parameter fiju/i, Ph = 2, /3„ = 4, Qh/CLm = 0.2. 
For Q M h = 0.275: f2 M = 0.69 h = 0.4 (solid line), 
Q M = 0.55 h = 0.5 (dashed line). For Q M h = 0.175: 
Q.M = 0.35 h = 0.5 (dotted line), « M = 0.26 h = 0.67 (dot 
dashed line) 

the neutrino rest mass, for one massive and two massless 
neutrino species is shown in Fig. 2. Increasing £Ih from 
0.1 to 0.3 decreases the power at galaxy scales by a fac- 
tor of about ~ 4 , at galaxy cluster scales by <~ 1.4 . 
The inverse behavior of these spectra near the maximum 



at k ~ 0.04 — O.ObhMpc^ 1 is caused by the existence of 
three characteristic timescales for the hot component: the 
epoch when it becomes non-relativistic, the period of colli- 
sionlcss damping and when the energy density of massless 
neutrinos equals the massive one. 

The power spectra for different f^A and fixed values 
for the parameters (3h , Pv,$Ih /^m are shown in Fig. 3. 
The horizontal shift of the spectra is caused by the later 
epoch of equality of the energy densities for the v and H 
particles. In different A models with fixed VIr the epoch of 
equality is different, in particular increasing the value of 

then Oc* decreases and z eq decreases too. The vertical 
shift is caused by the coefficient Ks/Qm in Eq.(10) and the 
amplification of CMB anisotropies in non-zero A models, 
which is included in the normalization parameter Shor- 
The change of £Ih in non-zero Oa models changes the 
spectra in a way similar to the f^A = case. 

The following property of spectra is suggested from our 
calculations: the models with the same component species 
{0h> fivi Qh/&m) but different Qm and h have the same 
form and localization of the maximum of P(k) when J7m^ 
is constant. This is shown in Fig. 4 for two values of I^m^: 
0.176 and 0.275. 

4. Large scale structure: predictions versus obser- 
vations 

4-1. Power spectra data and models 

In order to test our models we compare the linearly 
evolved P(k) with the one obtained from clustering data. 
Our power spectra P(k) have been computed for MLM 
models with a zero-baryonic content, in order to prop- 
erly compare with data we then correct the shape of 
the transfer functions according to a baryonic fraction 
tt b = 0.015//i 2 (see sect.4.4). 

We will make use of the dimensionless spectrum 
A 2 (fc) = k 3 P(k)/2-K 2 , calculated at z = 0, and com- 
pare it with the real space power spectrum of optical and 
IRAS galaxies obtained by Peacock (1997) and Peacock 
& Dodds (1994), taking into account redshift distortions 
and bias. Because we have performed a linear calculation 
we restrict ourselves to comparing A 2 (k) only in its lin- 
ear part fc^ Q.2hMpc~ l . For a linearly evolved spectrum 
the deviation from non-linearity because of gravitational 
clustering becomes significant at k£ O.ShMpc^ 1 . 

In Figs. 5 & 6 we show A 2 (k) for a set of MLM models 
with different values of Oa and Qh/Qm = 0.2 (Fig. 5), 
Qh/Qm —0.1 (Fig. 6). We assume one species of massive 
neutrinos. 

For each case we show two subcases: h = 0.5 and 
h = 0.7. For h = 0.7 we plot also the A 2 (k) of the ACDM 
model with ft a = 0.7, normalized to an rms mass fluc- 
tuation a(R = 8/i _1 Mpc) = cr 8 = 1 (see Eq.(19)). For 
this model we will use the notation A A (fc) for A 2 (fc). The 
observational data are shown as black circles and corre- 
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Fig. 5. a A (fc) is shown for MLM models with h = 0.5, (5h = 2, /3„ = 4, Qh/Qm = 0.2 and different values of the parameter: 
Qa = 0. (solid line), Qa = 0.31 (dotted line), ^a = 0.45 (short dashed line), Qa = 0.65 (long dashed line), Qa = 0.74 (dot-short 
dashed line). Black circles are the real space APM power spectrum obtained by Peacock (1997). b The same as in panel a, but 
for h = 0.7. Open squares represent the standard ACDM model with 57a = 0.7, normalized to ug = 1 



spond to the real space APM power spectrum for galaxy 
data (Peacock 1997). We make use of these data, rather 
than those of Peacock & Dodds (1994), because of the 
improved treatment for clustering evolution performed in 
the former paper. 

In order to constrain the parameter space for differ- 
ent models we have preferred not to do a \ 2 analysis be- 
cause of the correlation between different bins. Also be- 
cause of the statistical uncertainties we consider it more 
appropriate to make a qualitative comparison. We for- 
mally define a given model to be consistent with clustering 
data if the computed A 2 (A:) are within the 2a errors for 
fc< O.lhMpc- 1 . 

From Fig. 5a the model predictions for A 2 (fc) in the 
case of large values of fl\ (~ 0.3) are below the observa- 
tional data in the linear k region. 

This limit on J7a can be reduced if one considers h = 
0.7 (Fig. 5b). In this case the inconsistency in the linear 
regime is for J7a ^ 0.65. This shift in the spectrum follows 
from the dependencies discussed in the past section due 
to the increase in h. 

It is important to note the drastic reduction of small- 
scale power in MLM models because of massive neutrinos. 
In the parameter space the closest model to ACDM, for 
the set of MLM models of Fig. 5b, is Q A = 0.65. At high 
wave-number the spectrum of this model is well below the 
A A (fc) of the ACDM. For ease of comparison the ACDM 
has been normalized to as = 1, while the COBE normal- 
ization yields as — 1.1. On the other hand, for the MLM 
models considered, we always found <7s~ 1 (see sect. 4. 5). 



We thus conclude that MLM models have the property of 
removing the unpleasant feature of galaxy anti-biasing, in- 
voked for ACDM flKlypin, Primack fc Holtzman 1996| ) in 



order to fit clustering data. 

As we can see from Fig. 5 the spectra of MLM 
models with VLh/Qm = 0.2 do not contradict on lin- 
ear scales the spectrum reconstructed from observa- 
tions for fl M > 0.7(0.35) and h = 0.5(0.7). Fig. 6 
refers to Qh/Qm = 0.1 and consistency is achieved for 
n M > 0.55(0.25) and h = 0.5(0.7). 

This range of limits can be further constrained if one 
considers that non-linearity effects due to gravitational 
clustering will enhance A 2 (fc) well above our linear A 2 (fc) 
, for k<L QJlhMpc -1 (Peacock 1997). We then require that 
our linear spectra should be at least below the observed 
reconstructed spectrum in the high-fc region. Furthermore, 
for h — 0.7, we already know that the ACDM over pro- 
duces small-scale power, then the A 2 (k) of our MLM mod- 
els should not exceed the A A (fe) at high wave-numbers. 

The case ft = 0.5 does not give useful constraints, while 
for h = 0.7 we obtain 0.35< n M < 0.55 for Q h /^m = 0.2; 
0.25;$ tt M < 0.35 for Q h /&m =0.1. The constraints that 
we obtain are for a constant linear bias, with b ga i = 1. 

If the optical data are a biased tracer of the dark mat- 
ter then these limits should be changed according to the 
value of bg a i . However no simple scaling is possible for the 
obtained constraints because the value of b ga i will change 
for any model according to the value of S]« and Qh/Qm- 
In sect. 4. 5 we compare the constraints on Q,m that we 
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Fig. 6. a A 2 (fc) are shown for the same set of MLM models of Fig. 5a, but with a ratio Qh/^m 
a, but for h = 0.7 



0.1. b The same as in panel 



obtain from cluster abundances with those given by clus- 
tering data , according to the value of os ■ 

We do not attach a particular statistical significance to 
these lower limits on Qm, nevertheless in some cases we 
found these limits to be consistent with the constraints on 
Qm obtained independently from cluster abundances (see 
sect.4.5). 

4-2. Bulk motions 

Another constraint on dark matter models comes from 
the study of galaxy bulk flows in spheres around 



our position (see e.g., 


Kofman, Gnedin & Bahcall 1993J 


Stompor, Gorsky & Banday 1995 


Liddle et al. 1996b). 



Bcrtschinger et al. (1990) and Dekel (1994) give the aver- 
age peculiar velocities within spheres of radius between 10 
to 60 h Mpc after previously smoothing raw data with 
a Gaussian filter of radius Rf — 12h~ 1 Mpc. When the 
power spectrum is known then the rms peculiar velocity 
of galaxies in a sphere of radius R, corresponding to these 
data, can be calculated with the following expression 



confidence level of the rms predictions. A similar conclu- 
sion holds if ft, = 0.7. But taking into account the large 
error bars we must admit that these data do not rule out 
any of the models analysed here. We only conclude that 
models with CIm > 0.3 and ^h/^m < 0.2 are preferred. 

4-3. Cluster- cluster correlations 

The space distribution of rich clusters is a powerful tool 
for probing the power spectrum at intermediate and large 
scales. The first attempts to measure the two-point spatial 



autocorrelation function £ cc (r) ( Bahcall fc Soncira 1983 r 
Klypin fc Kopylov 1983| ) have shown that they are more 



V 2 {R) 



P(k)exp(-k 2 R 2 ' 



t , f )W 2 {kR)dk, 
(11) 

where W(kR) is the top-hat window function. 

The calculated predictions for the rms bulk motions for 
different MLM models with h = 0.5 are shown in Fig. 7. 
As we can see, decreasing Hm effectively reduces the bulk 
motion and the explanation of the observed data in low 
VLm models is problematic. So, in models with VLm < 0.3 
and flff/flM > 0.2 the observable data are above the 95% 



clustered than galaxies and ^(r) ss (ro/r) 1 with 7 = 
1.6 - 2.0 and r = 16 - 25 h^Mpc. 

The later analysis of other authors (see for example 
Postman et al. 1992; Olivier et al. 1993; Jing & Valdarnini 
1993) has confirmed that £ cc is well fitted by the same 
expression with 7 = 1.8 and ro w 20 h Mpc. The impor- 
tant conclusion, which follows from numerous studies of 
this problem, is the existence of a positive long distance 
correlation of rich clusters of galaxies out to 50 h~ 1 Mpc. 
For a Gaussian random density fluctuation field the corre- 
lation function of peaks at large separations is calculated 



with the following equation (Bardeen et al. 1986): 



Mr) 



f% / P(k)k 2 W 2 (kR c 
2tt 2 Jo 



sin kr 



kr 



-dk, 



(12) 



where W{kR c ) is the window function, which filters out in 
the density field the structures on scales larger than R c , b c 
is their biasing parameter, which takes into account the 
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Fig. 7. The bulk motion V rm a(R) is shown for MLM models 
with h = 0.5, fa = 2,p u = 4, Q H /Sl M = 0.2 and different 
values of the Qa parameter: Qa = (solid line), Qa = 0.31 
(dotted line), Qa = 0.45 (short dashed line), Qa = 0.65 (long 
dashed line), Qa ~ 0.74 (dot-short dashed line). The solid 
circles correspond to Dekel(1994), open squares to Bertschinger 
et al.(1990), solid squares to Courteau et al. (1993) 




0.01 r 
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100 
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Fig. 8. The cluster-cluster correlation functions £ cc are shown 
for MLM models with h = 0.5, f3 H = 2,/3 v = 4, Q h /Qm = 0.1 
and different values of the Qa parameter: Qa = (solid line), 
Qa = 0.31 (dotted line), Qa = 0.45 (short dashed line), 
Qa = 0.65 (long dashed line), Qa = 0.74 (dot-short dashed 
line) 



statistical correlation of peaks above a given threshold. 
The biasing parameter is defined by the expression 

b c =<D> /<i(R c ) + 1, (13) 

where the effective threshold level < v > is given by 

70 



< v >- 



dv 



l-Y 



t{v/v t )N pk (v)/n c (14) 



with 7, 9, the threshold function t(i>/v t ) and the dif- 
ferential number density N p k(is) are defined according 
to equations (4.6a), (6.14), (4.13) and (4.3) of Bardeen 
et al. (1986). The Gaussian filter radius correspond- 
ing to the mass (M > 5 10 14 /i _1 M Q ) of a rich clus- 
ter is R c = (M/4.35 lCP^Af^ 1 ) 173 /*" 1 Mpc. The 
observed number density of Abell clusters with rich- 
ness R > 1, nf s = (5.7 ± 0.5 ) 10^ 6 h 3 Mpc~ 3 (see 
Zamorani et al. 1991 |Bahcall 1988 for a review), is used 



for the determination of the peak height ut 



dvt(v / v t )N pk {y) = n° 



bs 



(15) 




0.1 - 



0.01 r 



1E-3 



100 



The mean height of peaks for which the rich clusters 
of galaxies are formed is then 



< v >= 



dvt(v/v t )vN pk (u)/n c 



(16) 



r (h Mpc) 

Fig. 9. The same as in Fig. 8 , but for Qh/Qm = 0.2 
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The cluster-cluster correlation function calculated in 
this way for the standard CDM scenario strongly con- 
flicts with the one observed , because it becomes nega- 
tive after ~ 40 h~ 1 Mpc. For the standard MDM mod- 
els with Hh < 0.3 the situation is better (£ cc become 
negative at ~ 50 h~ 1 Mpc) but positive correlations at 
> 50 h~ 1 Mpc are not explained by them either (see for 



example Holtzman fc Primack 1995; Novosyadlyj 1994) 



The positive correlations at > 50 h 1 Mpc would be 
explained in the cosmological scenario with Qm = 1 



large scale power ( 


Bardccn, Bond & Efstathiou 1987 , 


Novosyadlyj & Gnatyk 1994 




Novosyadlyj 1996|). But 



such spectra are now ruled out by the data for AT/T at 



degree angular scales ( |Schustcr et al. 1993|) because their 



prediction for the rms AT/T is higher than the 95% c.l. 
of the experimental upper limit. 

An attractive possibility for avoiding this problem 
is to consider MLM models normalised to the COBE 
quadrupole. The correlation functions of rich clusters of 
galaxies calculated as described above agree with obser- 
vational data and are positive far beyond 50 h^Mpc 
(Figs. 8 & 9). As we can see, for MLM models with 
0.3 < n M < 0.7, n H /n M = O.l ~ 0.2 and h = 0.5 the pre- 
dicted autocorrelation functions of rich clusters are within 
the limits of the observed error bars at large distances. A 
comparison for other h gives the following constraints for 
fl M h: 0.13 < n M h < 0.35. 

The biasing parameters of rich clusters of galax- 
ies for all models are in the range 3.3 — 4.6 and 
are in the same range as the values obtained with 



different methods from observations ( 


Lyndcn-Bcll 1991, 


Plionis & Valdarnini 1991 




Plionis 1995 


) . The minimal 



Q]\jh is constrained also by the moment of turn around of 
density fluctuation peaks which are associated with rich 
clusters. Indeed, the requirement that the total cluster 
mass collapsed before z = (not only the central region 
but also the frontier areas ) requires < v > <j c > 1.06 
(Bardeen et al. 1986), which is satisfied for models with 
flujh > 0.15 when Ojj/Ojvf = 0.1 and D,Mh > 0.18 when 

n H /n M = o.2. 



4-4- Cluster mass function 

The observed abundances of clusters of galaxies is a pow- 
erful discriminant for different models of dark matter. We 

(1974, PS) formula 



will use the Press & Schechter 1974 



to compute the number density 



JV(> M,z) 



n(M',z)dM', 



(17) 



of virialized objects with mass greater than M. Accord- 
ing to PS, the comoving number density n(M, z) of halo 
masses in the interval M, M + dM , is analytically related 



to the power spectrum by 



n(M,z) 



Sc. 



7T M 2 a{M,z) 



dlna 



dlnM 



-S 2 c /2a 2 



(18) 



here Pb{z) is the background density, a(M,z) is the rms 
mass density fluctuation and S c is the threshold parame- 
ter. The rms a value is defined as 

f°° rlk 

a 2 (M,z) = D 2 (z) — A 2 (k)W 2 (k,R), (19) 
Jo k 

where W is a window function and D is the linear growth 
factor, -D(O) = 1. The relation between M and R de- 
pends on the choice of W, for a top-hat function R = 

According to linear theory, 5 C — 1.686, for a top- 
hat window in an Qm = 1 universe. For the more 
general case VLm < 1, 5 C can be derived analyt- 
i cally and it has a ve ry weak dependence on Qm 



( |Eke, Cole fc Frenk 1996[ ). On cluster scales, density fluc- 
tuations are well described by linear theory and so Eq.(18) 
is thought to be a good approximation to the true num- 
ber density. N-body simulations have been used by var- 



ious authors ( 


Efstathiou et al. 1988; 


Lacey & Cole 1994 


Eke, Cole & Frenk 1996 ) for an extensive check of the va- 



lidity of the PS predictions. In particular Eke et al. (1996) 
found that Eq.(18) agrees well with N-body results for 
CDM models in a flat universe. 

For a top-hat choice, the best-fit to N-body results is 
5 C = 1.7 ± 0.1, while for a Gaussian window the 8 C thresh- 
old is more sensitive to the shape of the power spectra 



(Lacey & Cole 1994 ). In what follows we will take a top- 
hat window function W and S c = 1.686. Our results will 
be compared with those of Ma (1996) , who has made a 
similar choice. Cluster abundances have been computed 
for MDM by Ma (1996) ; Bartlett & Silk (1993); Liddle et 
al. (1996c); Bahcall, Fan & Cen (1997). Clusters of galax- 
ies are rare objects and therefore are very sensitive to the 
value of a . Our models have two limiting cases: when 
fl A -> ( MDM models) or when fl H /^M -* , which 
corresponds to ACDM models. For these limiting cases we 
have compared our integrations with published analytical 
spectra and found that the greater discrepancies are at 
high wavenumbers fc<; lhMpc^ 1 and are ~ 10% of the 
analytical P(k). 

These differences are mainly due to the absence of 
baryons in our computations. For the length scale of in- 
terest to us, R ~ 8 — lQh~ 1 Mpc and the resulting a's will 
have similar relative errors. Because of the exponential in 
Eq.(18), the cluster number density is strongly affected by 
these errors. We have therefore decided to include baryons 
in our computations of linear spectra. In fact we have ob- 
tained the new spectra by taking the computed transfer 
functions T(k) and interpolating them linearly over a grid 
of values q = kj^uh 2 . 
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Fig. 10. The cluster abundance N(> M) is shown at z = 
for MDM models with Qh/^m = 0.2, h = 0.5 and different 
values of the Ph and f3 v parameters. The values chosen are: 
Ph = 2 f3 u = 4 (solid line), p H = 6 p v = (dotted short line), 
Ph = 4 p v — 2 (short-dashed line). Data points have been 
taken from Ma (1996, see text). Black squares correspond to 
the integration of Ma for the model Ph — 2 p v = 4. The linear 
spectra used for the calculation of N(> M) are normalized 
according to the COBE 4-yr data 



For the same wavenumber k as in the original com- 
putation, the new T{k) is obtained according to the pre- 
scription of Sugiyama (1995) : q — > q/exp(—Qb — Q&/Qm)- 
We have taken = 0.015 /h 2 . This procedure works well 
for CDM transfer functions because, after recombination, 
baryons are caught by the CDM component and their per- 
turbations grow together. For spectra with a hot compo- 
nent the ratio Th/Tcdm to a first approximation does 
not depend on Qm and the prescription can be applied to 
the total transfer function. For MDM Liddle et al. (1996c) 
have applied a similar procedure to the transfer functions 
of Pogosyan & Starobinsky (1995), who do not include 
baryons in their calculations. They found that the appli- 
cability of the procedure requires 0.1 and hft 0.5, a 
range of limits that we never consider. 

In order to compare our results with observations we 
take the number density of clusters at two mass ranges 
from the work of Ma (1996). These data points correspond 
to clust ers with X-ray tempera ture greater than 3.7 and 
7 keV ( Henry fc Arnaud 1991 ) . For the first point White 



mass of M = 4.2 10 14 M Q assuming an isothermal model. 
The second point is taken from Liddle et al. (1996d) who 
have used the N-body hydro simulations of White et al. 
(1993b) and found M = (1.23 ± 0.3)W 15 h- 1 M Q for the 
virial mass of a cluster with X-ray temperature kT = 7 
keV in a critical universe. 



o.ooi 



o.oooi ? 




10" 



10 15 10" 
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et al.( 1993a) estimated the upper limit in mass from the 
cluster velocity dispersion, while for the lower limit the 
X-ray temperature of 3.6 KeV has been converted into a 



Fig. 11. The same as in Fig. 10 but for MDM models with 
h = 0.5, p H = 2,P„ = 4, and Q h /^m =0.1 (solid line), 
Qh/Qm = 0.2 (dotted line), Qh/Qm = 0.3 (short dashed line) 



In Fig. 10 we plot the present mass functions of clusters 
for SMDM models with n H = 0.2, h = 0.5, ft b = 0.06. The 
figure shows the dependence of N(> M) on the number of 
species of massless and massive collisionless particles. We 
assume that the total number of species of both particles 
is equal to 3. We can see a weak dependence on these 
parameters. The mass function in SMDM models with one 
species of massive neutrinos is above that for models with 
three massive neutrinos. Keeping Qh fixed, and decreasing 
the neutrino mass, then N(> M) decreases too. It can be 
seen that the model with Ph = 4, advocated by Primack 
et al. (1995), fares much better. 

Fig. 11 shows the mass function N(> M) for SMDM 
models with different values of Vl H = 0.1, 0.2, 0.3 and fixed 
values of h — 0.5, Ph — 2, /3„ = 4. The black squares corre- 
spond to the £Ih = 0.2 case treated by Ma (1996, figure 8 
top right). With an increase of fin , the mass function de- 
creases. We can observe also small differences in the shape 
of the functions. 

As we can see, for the considered range of parameters 
SMDM models do not fit cluster abundances, if we adopt 
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the COBE normalization and spectral index n = 1 with 
no gravitational waves. The fin = 0.3 case is marginally 
consistent, and it was the one considered by Bartlett & 
Silk (1993). An improvement in the fit can be obtained 
with the introduction of a small tilt in the initial spec- 
trum ( n ~ 0.9 ) and/or a tensor contribution to COBE 
anisotropies. This point has been discussed in detail by 
Ma (1996) and we do not consider it here. 

For MLM models the data points in Fig. 11 should be 
calculated according to the model, because the conversion 
from temperatures to masses depends on £Im- For this 
reason we next consider the cumulative cluster temper- 
ature function predicted for different MLM models, and 
compare it with observations. 

4-5. Cluster temperature function 

Identification of clusters of galaxies in the optical band 
is subject to the problems of foreground/background con- 
tamination. Projection effects can also undermine mass 
determination through virial analysis ( Frcnk ct al. 199lj[ 
Dckcl ct al. 1989] ). On the other hand c lusters are also 
strong X-ray sources ( McKee et al. 1980 ), their emission 
does not suffer from these problems and clusters can be 
reliably identified. 

The X-ray emission of galaxy clusters has two phys- 
ical observable- related quantities: the luminosity and 
the temperature. During cluster collapse, the gas is 
shock heated to the virial temperature, then it ap- 
proaches an isothermal distribution in virial equilib- 
rium. For the gas temperature T g one should have 
T g cx aM 2 / 3 . This relation has been confirmed by nu- 



merical simulations ( Evrard, Metzler fc Navarro 1996| ; 



Navarro, Frenk & White 1995) with a very small disper- 
sion for the coefficient a and different DM models. 

Then the T — M relation allows us to connect the PS 
equation (18) to the cluster X-ray temperature function 
(XTF). We take the Eke, Cole & Frenk (1996, hereafter 
ECF) relation for an isothermal gas: 



hTfj, 



7.75 ( 6.. 



\5X 



M 



(1+*) 



Ac 
178 



keV, 



(20) 



here (3 is the ratio of the galaxy kinetic energy to the 
gas thermal energy, X is the hydrogen mass fraction, A c 
is the ratio of the mean halo density within a virial radius 
to the critical density at the corresponding redshift. We 
assume X = 0.76, (3 = 1 and A c = 178 (ECF). For Q M < 
1, A c can be derived analytically and is well approximated 
by A c = 178ft°/ 5 . 

For a specified model Eqs. (18-20) allow us to com- 
pute the cluster XTF N{> kT). Then this function 
can be compared with present data and can be used 
to constrain different DM models ( Bartlett fc Silk 1993 
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Fig. 12. We show the temperature function of clusters 
N{> kT) for MLM models with h = 0.5, f3 H = 2,/3„ = 4, 
Qh/£Im = 0.2 and different values of the Qa parameter: 
Ov = 0. (solid line), Ov = 0.31 (dotted line), fl A = 0.45 
(short dashed line), Oa = 0.65 (long dashed line), S1a = 0.74 
(dot-short dashed line). The histogram is obtained from the 
X-ray data set of Henry and Arnaud (1991). Open squares are 
three selected temperatures for which error bars have been es- 
timated with a bootstrap procedure. Black circles refer to the 
mass function of Bahcall and Cen (1993). The horizontal bar 
is from White et al.( 1993a) 



ECF; Viana fc Liddlc 1996 ). Current observations from 
Einstein and EXOSAT satellites have been used to de- 



termine the loc al cluster XTF ( Henry fc Arnaud 1991 
Edge ct al. 199CD - Our study of the cluster XTF will make 



use of the Henry & Arnaud (1991) data. In particular we 
will closely follow the analysis of ECF, this will make it 
easy to compare of our MLM predictions with previous 
result for ACDM models (ECF). 

The estimated cumulative cluster XTF is computed 
according to 



N(> kT) 



Ti>T 



i/v„. 



(21) 



here V ma x,i is the maximum volume at which the i — th 
cluster can be detected for a specified flux limit (Fx = 
3 10~ 11 ergsec~ 1 cm~ 2 ) in the 2—10 Kev range. 

The cumulative XTF obtained in this way is shown as 
the solid line histogram in Fig. 12. Error bars have been 
found with a bootstrap procedure, applied to the origi- 
nal sample of 25 clusters, at three different temperature 
bins Tj. The resulting ler amplitudes 5j are plotted as 
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open squares in Fig. 12, at the three corresponding tem- 
peratures. The Sj arc defined as <5| =< (log 10 N boot (> 
kTj) - log 10 N boot {> kTj)) 2 >, the average being over the 
bootstrap ensemble. 

The estimated A^(> kT) agrees with the one obtained 
by ECF (see their Fig. 3), the procedure being the same, 
and is consistent with the differential one of Henry & Ar- 
naud (1991), although ECF pointed out some errors in the 
original estimate that almost cancel each other. 

The horizontal bar on the left upper part is the first 
data point of the previous iV(> M) figures. The point 
was obtained from White et al. (1993a) precisely from the 
temperature function of Henry & Arnaud (1991) and it 
proves the consistency of our procedure. The small offset 
along the kT axis can be attributed to the use of the 
differential or cumulative temperature functions (ECF). 

A different test is to compare our N(> kT) with the 
cluster mass function of Bahcall & Cen (1993). Their best 
fit N(> M) is plotted as black circles for three different 
temperatures. In order to convert masses into tempera- 
tures we have used Eq.(5) of Bahcall & Cen (1993). Error 
bars represent the la dispersion of their Eq.(5). As can 
be seen, there is a substantial agreement between the two 
estimates. 

In comparing the predictions of our models with the 
observed XTF one must be aware of possible uncertain- 
ties in the mass-temperature relation. Because of the steep 
decline of the XTF with temperature, even a small disper- 
sion can have drastic effects. In Eq.(20) the parameter j3 
represents the ratio between the virial and gas tempera- 
tures (we have assumed an isotropic isothermal profile). 
The (3 parameter can be calculated either from the data 
or through numerical simulations. 

Current available data are consistent with f3 ~ 1 
(|Edge fc Stewart 199l| |Squires at al. 1996|). The 



dynamical simulations of Navarro, Frenk & White (1995) 
give a result of /3 ~ 1 and have been used by ECF , who as- 
sume (3 = 1±0.1. This scatter around unity for j3 ( ~ 10%) 
has been confirmed in more recent work by Eke, Navarro 
& Frenk (1997, see also Evrard, Metzler & Navarro 1996), 



who have performed a set of N-body hydrodynamical sim- 
ulations to investigate the X-ray evolution for a set of 
clusters in a low-density flat CDM cosmology. 

Thus we have reliably assumed /3 = 1 in order to con- 
vert masses into temperatures. The theoretical XTF has 
then been computed from Eq.(18), for MLM models with 
different J7a and a fixed ratio VLh/Qm = 0.2. We nor- 
malize the final spectra according to the COBE data (see 
sect. 2. 2 ). The results are shown in Fig. 12. 

The inclusion of a A term clearly alleviates the prob- 
lem for MDM and brings the models in good agreement 
with the data. The solid line is the limiting case Q\ = 0, 
which corresponds to the standard MDM. In order to fit 
the present cluster abundance, this model would require 
a value j3 — 2, clearly inconsistent with present estimates. 
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Fig. 13. The same as in Fig.12, but for Q h /^m = 0.1 



From Fig. 12 the best range of fl\ for fitting the data 
is fl\ — 0.3 — 0.5. Fig. 13 shows the same plot of Fig. 12 , 
but for CIjj/CIm = 0.1. In this case the best range for Q\ 
is ~ 0.4 — 0.6. In order to constrain a particular model we 
can now compute the x 2 quantity 



(22) 



3 = 1 k=l 



hydro- where j is the bin index, yj — (log 10 (Nj) 



log 10 (iVf)), 

yj — N(> kTj) and N* h is for a particular model from 
Eq.(18). The covariance matrix Cjk =< SjSk > takes 
into account the correlations between different temper- 
ature bins. According to ECF these correlations are not 
negligible, but the models which minimize \ 2 do not de- 
pend strongly on their inclusion in Eq.(22). Thus ECF 
take for Cki its diagonal form. However the models that 
we consider have different spectra and we have chosen to 
keep the whole matrix for the minimization of Eq.(22). 

From Eqs. (18-19) \ 2 is a function of the power spec- 
trum constant A or, cquivalcntly, of the rms mass fluc- 
tuation erg = a(R — 8/j -1 ). From the minimization of 
Eq.(22) the formal error on a% is 5 %, but uncertainties 
on the other parameters will affect a% too. From N-body 
integrations ECF assume a scatter of 4 % for 5 C . This is 
not surprising because for the length scale of interest to us 
we can neglect tidal forces and assume that the collapse 
is spherically symmetric. Other sources of errors are the 
scatter in j3 ( ^ 10 % ) , the sam ple completeness ( 90 % 
at Fx = 3 10~ 11 ergsec~ 1 cm~ 2 , Lahav et al. 1989] ) and 
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Fig. 14. a as is shown as a function of Om for MLM models with h = 0.5, (3h = 2,/3„ = 4, Qh/Qm ~ 0.2. The continuous 
line is the best-fit to the estimated cluster abundance, the dashed line is from the COBE 4-yr data. The thin lines show the 
statistical uncertainties, quoted in the text, b The same as in panel a, but for Qh/^m = 0.1 



errors in the measurement of temperatures. For the latter 
error, ECF analysis gives an upper limit of 1 %. Sum- 
ming all these errors in quadrature the final dispersion for 
as is about 10 %, twice the statistical error. 

For a given model we can now estimate, from Eq.(22), 
the erg which is consistent with the estimated cluster abun- 
dance, and compare it with the as obtained from COBE 
data. These two values for as will in general not coincide, 
however there will be a set of values of cosmological pa- 
rameters for which they will be in the same range. The 
estimated uncertainties allow us to judge the reliability 
of the overlap. For the COBE normalization we ha ve as- 
sumed a 7 % statistical error ( Bunn fc White 1997 ). 

We show in Fig. 14a the as ( continuous line ) which 
is obtained from the minimization of the x 2 quantity. The 
dashed line is from COBE. The two sigmas are plotted 
as function of the cosmological parameter £Im- The set 
of models is for h = 0.5 and VLh/^m — 0.2. Thin lines 
represent the assumed uncertainties. 

From Fig. 14a the best range of £Im which fits the data 
is 0.55< VL M < 0.75. The standard MDM model (fi M = 1) 
is rejected at the 2cr level. Fig. 14b shows the same plot 
as in panel (a), but for H,h/^m = 0.1. In this case the 
best range for Om is 0.45^5 £1m~ 0.65. In Fig. 15 we con- 
sider the same set of models but for h = 0.7. We ob- 
tain 0.3< fi M < 0.5 for n H /n M = 0.2 and 0.3< fi M < 0.4 
when £Ih/^m = 0.1. 

According to the value of h we can then summarize the 
following constraints for MLM models with one species 
of massive neutrinos and VLh/Qm < 0.2: 0.5^$ f^M~ 0.8 
(h = 0.5) or 0.3< Q M < 0.5 (h = 0.7). The correspond- 



ing allowed values of as are : a s — 0.65 ±0.1 (h = 0.5) 
or as — 0.8 ± 0.1 (h = 0.7). A comparison of these con- 
straints with the clustering data of sect. 4.1 must take into 
account the biasing factor. A rescaling of A 2 (fc) according 
to bgai — I/eg shows that consistency with cluster abun- 
dance is marginal for models with VLh/SIm = 0.2. For 
h = 0.5 and n H /fl M = 0.1 we obtain 0.55< n M < 0.65 , 
while for h = 0.7 the limits on SIm are 0.35^ £Im^ 0.4. 

The constraints on 17m that we obtain follow from 
comparing our results with the present cluster XTF. We 
have not considered possible constraints that might be 
given by considering cluster evolution. The evolution of 
cluster number density with redshift is a powerful tool 
for discriminating among different cosmologies. This fol- 
lows directly from the different growth of density fluctua- 
tions in different cosmologies. The fluctuations in models 
with CIm < 1 grow very slowly and structures will experi- 
ence little evolution at recent times. On the contrary, for 
= 1, the higher growth rate of density fluctuations 
implies for galaxy clusters a much stronger evolution at 
late redshifts. 

Bahcall, Fan & Cen (1997) compared the results for 
cluster abundance, using large-scale N-body simulations, 
with recent data at z ~ 0.5 - 1. They found that MDM 
models are ruled out at the 2a level, while a ACDM model 
with £Im = 0.34 ± 0.13 h = 0.65 can consistently fit the 
data. These simulations show that the cluster evolution 
rate is strongly model dependent. We accordingly do not 
attempt here to perform a linear analysis, applying Press- 
Schechter theory, to compute the redshift evolution of the 
cluster number density for MLM models. 
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Fig. 15. a as is plotted versus Qm for MLM models with h — 0.7, (3h — 2, /3„ — 4, Qh/^m = 0.2. The meaning of the lines is 
the same as in Fig. 14. b The same as in panel a, but for JJg/fljf = 0.1 



We suggest then that large scale numerical simula- 
tions can be used to compare the cluster evolution for 
MLM models with recent data. These tests are likely to 
strengthen or falsify the models, given the already narrow 
window of allowed values for the cosmological parameters 
arising from the linear tests applied here. 

4-6. Damped Lyman-a systems 

An important test for dark matter models with a massive 
neutrino component is given from the observations of ob- 
jects at high redshifts. Due to free-streaming effects the 
dark matter power spectrum will be severely damped on 
small scales, thus making formation of early objects more 
difficult than in a model without the hot component. The 
most important class of such objects are damped Lyman-a 
absorption systems. 

These objects have a high column density of neutral 
hydrogen (Nhi<^ 10 20 cm~ 2 ) and are detected by means 
of absorption lines in quasar spectra (Wolfe 1993). Ob- 
servations at high redshift have lead to estimates of the 
abundance of neutral hydrogen in damped Lyman-a sys- 
tems (Lanzetta, Wolfe & Turnshek 1995; Storrie-Lombardi 
et al. 1996). The latter authors have analyzed spectra at 
z = 3 and z — 4 , we will make use of the z — 4 data, which 
gives the strongest constraint on the primordial spectrum. 
According to Storrie-Lombardi et al. (1996), the cosmo- 
logical mass density of neutral hydrogen gas Q g is given 

by 



n g (z = 4) = (o.oon ± o.ooo2)/r 1 < 



n 



n M {z = *y 



(23) 



with the original data being given in units of the criti- 
cal density, the square root taking into account different 
cosmologies. 

The standard view is that damped Lyman-a systems 
are a population of protogalactic disks (Wolfe 1993), with 
a minimum mass of M > lO lo /i _1 Af0 (Hachnclt 1995). 
The fractional density of collapsed objects of minimum 
mass Molas is then 

n DLAS (z = A) = " g(z = 4) =(0.069 ±0.021) -?L 

Jgas^ t'b JgcLS 



9. 



n M (z = 4) 



(24) 



where f gas is the fraction of neutral hydrogen and 
Q, b = 0.016/i~ 2 flCopi, Schramm fc Turner 19950 . In Eq. 



(24) the error in fif, has been added in quadrature (Liddle 
et al. 1996c). A conservative assumption is f gas = 1, but 



recent hydrodynamical simulations (Ma et al. 1997) have 



claimed f gas < 0.1. We will consider constraints on our 
models arising from both of the limits on f gas . 

The theoretical counterpart of Eq.(24) can be found 
using the Press-Schechter (1974) theory and is given by 



n 



DLAS 



(> M, z) = crfc 



V2a(M, z) 



(25) 



where crfc is the complementary error function and 
a(M,z) is defined according to Eq.(19) using a top-hat 
window function. Particular care must be taken when 
inserting the minimum mass into Eq.(19) because for 
damped Lyman-a systems M = Mdlas is well below 
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From Fig. 16 one can infer that for h = 0.5 MLM mod- 
els with CIh/^m = 0.1 are not strongly constrained by the 
available data for damped Lyman-a systems. For h = 0.5 
and Qr/^m — 0.2(0.1) we obtain the following lower lim- 
its : 0.8(0.4) < fi M . For n H /Q M = 0.2 this lower limit on 
fijif do not overlap with the upper limit given by cluster 
abundances and the model is inconsistent. For h = 0.7 we 
obtain from Fig. 17 0.4(0.2)< M if n H /Q M = 0.2(0.1). 
This lower limit is consistent with the upper limits from 
cluster abundances and other constraints for the same set 
of models. 

These constraints on Qm have been obtained for a 
minimum mass of M = 10 10 ft. _1 M Q , the role of possi- 
ble uncertainties on the limits for Qm can be considered 
by decreasing this lower limit by an order of magnitude. 
For M — 10 9 h~ 1 M ( ? ) we obtain new constraints on Am 
which do not differ in a relevant way from those shown in 
Figs. 16 & 17. 



Fig. 16. a(M,z = 4) is shown as a function of Qm for 
h = 0.5. The continuous line represents the 95 confidence level 
for the the lower limit obtained from observational data using 
Press-Schechter theory. The bottom line is for f gas = 1 and the 
top line refers to f gas = 0.1. The o(M,z = 4) for MLM mod- 
els are plotted as dashed lines for different values of the ratio 
Q,h/Qm- The long(short)-dashed line is for Qh/^Im = 0.2(0.1) 



the neutrino clustering scale at z = 4. We correct Mdlas 
according to M DLAS -» 10 10 /i" 1 (1-O h /^m) _1 M q (Lid- 
die at al. 1996c). The collapsed state of damped Lyman-a 
systems is uncertain and a conservative assumption is that 
these systems have collapsed along two axes. We accord- 
ingly take for the threshold parameter the value S c = 1.5 
(Monaco 1995). 

We conservatively assume the theoretical predictions 
(25) to be bounded from below by Eq.(24), for SI d las 
we take the 2a lower limit. The rms mass fluctuation 
<r(M, z = 4) can now be constrained from below using 
Eqs. (24-25). We show these limits as a function of VIm in 
Figs. 16 & 17, for h = 0.5 and h = 0.7, respectively. 
In these figures the bottom solid line is for f gas = 1, 
while the top solid line refers to f gas = 0.1. This value 
gives a much more severe constraint on a(M, z = 4) 
than for f gas = 1. The theoretical values of <r(M, z = 4) 
are shown for SIh/^-m = 0.1 ( short dashed line ) and 
Qh/^m = 0.2 ( long dashed line ). At a given Q M the 
value of a(M,z = 4) decreases as the ratio £L H /£l M in- 
creases. This behavior of cr(M, z = 4) follows because an 
increase of implies, for a given normalization and red- 
shift, a reduction in the available power at small scales. 
The effect becomes more pronounced for Qm — * 1 
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Fig. 17. The same as in Fig.16, but for h = 0.7 

We have considered f gas = 1 but if we take f gas = 0.1 
then none of the MLM models for h = 0.5 survives the con- 
straints from damped Lyman-a systems. The model with 
Qh/Qm = 0.2 is excluded and SIh/^m = 0.1 requires 
n M Z o.8, a minimum value excluded at the la level by 
cluster abundances. Also MLM models with h = 0.7 are 
severely constrained if f gas = 0.1. In this case we have 
0.8(0.4) for Qh/^m = 0.2(0.1). Thus the model 
with CIh/Qm = 0.2 is totally inconsistent with cluster 
abundances but Qh/^m = 0.1 is still within the la range. 
We conservatively take the lower limit on Qm arising from 
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fgas = 1. Numerical hydrodynamical simulations have 



tested only MDM models ( |Ma et al. 1997[ ), hydro simu- 
lations with MLM spectra are clearly required in order to 
obtain a tight limit on f gas - 

5. Conclusions 

We have discussed linear clustering evolution for a set of 
spatially flat MDM models with a cosmological constant. 
We have not considered the role of gravitational waves or 
of a possible tilt in the primordial spectrum. In order to 
restrict the range of allowed values for the cosmological 
parameters, we have applied linear perturbation theory to 
compare the predictions of our MLM models with a set 
of linear data. The models considered had one species of 
massive neutrinos and a scale-invariant spectrum. 

Linear calculations for the dimensionless power spectra 
A 2 (fc) have been compared with the reconstructed real- 
space power spectrum of Peacock & Dodds (1994), Pea- 
cock (1997). Because we have considered linear spectra 
evolution a comparison of clustering is meaningful only 
at low k (JS 0.2hMpc~ r ) A substantial agreement is ob- 
tained for those unbiased models with fln/^M^ 0.2, and 
n M > 0.6 (h = 0.5) or 0.3< ft M < 0.5 (h = 0.7). 

The computed linear power spectra have been filtered 
with a Gaussian window of radius Rf — 12h~ 1 Mpc to cal- 
culate the peculiar velocity field. A comparison of the rms 
bulk motion with POTENT data does not lead to strong 
constraints : f^A^S 0.7 and fln/^M^ 0.2. The cluster- 
cluster correlation functions calculated for MLM spectra 
in the framework of Gaussian random density fluctua- 
tion field explain the observed positive correlations at > 
50 h~ 1 Mpc and are within the limits of the error bars for 
that observed at large distances when 0.13 < Am^ < 0.35. 

Of the considered tests the most important turns out 
to be that based on cluster abundances, for which the 
present X-ray cluster temperature function can put strong 
limits on the initial spectrum amplitude, or conversely on 
a 8 . Because of the strong dependence of N c i(> kT) on 
temperature, the conversion of X-ray temperatures into 
masses is critical . From the discussion of Section 4.5 this 
conversion can be obtained with a dispersion of less than 
10 %. We have then used the Press & Schechter (1974) 
formalism to compute the cluster number density for dif- 
ferent MLM models. Consistency with present data can 
be achieved for COBE normalized models only for the 
following range of values (at la level): SIb/^m < 0.2, 
h = 0.5(0.7) and 0.45(0.3) < Vt M < 0.75(0.5). In corre- 
spondence with these limits, the allowed values of as are 
: cr 8 ~ 0.65 ± 0.1 {h = 0.5) or a s ~ 0.8 ± 0.1 {h = 0.7). 
These limits suggest that MLM models require a moder- 
ate amount of bias to be consistent with clustering data. 
For these range of values consistency with clustering data 
yields 0.55(0.35) < ft M < 0.65(0.4) for h = 0.5(0.7) and 
Qh/^M — 0.1. Models with Hh/^m — 0.2 are almost 
inconsistent. 



One can also compare these constraints with those 
given by considering cluster evolution. The evolution of 
galaxy clusters has recently been used also for determin- 
ing the cosmological parameters as and fijw. Fan, Bah- 
call & Cen (1997) found a s = 0.83 ± 0.15. This range of 
values is consistent with the as required for MLM with 
h = 0.7 (Fig. 15), and marginally consistent if h = 0.5. At 
the 2a level none of the MLM models considered here is 
ruled out. Finally Henry (1997) has used recent data for 
the evolution of the cluster X-ray temperature function 
to determine f2 m aud as ■ For a flat cosmology he obtains 
!]m = 0.55 ± 0.17 and as = 0.66t^i7j m cl° se agreement 
with our findings for MLM models. 

The observed abundance of damped Lyman-a systems 
has been used to obtain lower limits on the rms mass 
fluctuation a(M) at z = 4 using Press-Schechter theory. 
These limits yield the following constraints on Qm for 
MLM models with tt H /Q M < 0.2 : 0.8 - 0.4(0.3)< Q. M 
for h = 0.5 (h = 0.7). These values have been obtained 
assuming the fraction of neutral hydrogen to be unity. If 
this fraction is close to 10%, as suggested by numerical 
hydro simulations, then the constraints for MLM models 
become much more severe. In this case consistency with 
both damped Lyman-a and cluster abundances can be 
achieved only for h = 0.7 , Qm — 0.4 and SIh/^m = 0.1. 

In order to further restrict our models we can also con- 
sider present observational constraints on cosmological pa- 
rameters. The cosmological constant is constrained to be 
< 0.66 at the 95 % confidence limit by QSO lensing 
(Kochanek 1996). A tighter restriction comes from the re- 
cent work of Perlmutter et al. (1997) on SN la , who give 
fl A < 0.51 still at 95 % c.l. . 

Recent Hipparcos data ( Feast & Catchpole 1997, 
Rcid 1997) have brought down the Cepheid scale dis- 
tance, thus reducing the estimated Globular Cluster age 
to to > \2Gyr and ft to a mid term h = .65. These values 
still require a cosmological constant, but not as high as in 
ACDM models. A value of fl\ > 0.4 is needed to satisfy 
the new age constraint. 

An upper limit on Qm can be obtained from the esti- 
mated baryonic content of galaxy clusters. If clusters are 
massive enough to represent a fair sample of the total mat- 
ter co ntent in the universe, as numeric al simulations con- 
firm ( Evrard, Metzler fc Navarro 1996| ), then the baryon 
fraction f b = (0.06±0.003)fr~ 3 / 2 (Evrard 1997) should be 
close to its universal value. 

Thus the standard Big Bang Nucleosynthesis value 
(0.008 -0.024)ft" 2 (popi, Schramm fc Turner 1995[) 



can be used to infer Q,m = Sib/ h — (0.25 ± 0.15)h~ 
For h = 0.65 one obtains Qm < 0.5. For a flat model this 
limit already overlaps the lower limit from SN la. This is 
the main argument for a cosmological constant, possible 
counter arguments like a magnetic field pressure or den- 
sity inhomogeneities, which can lead to underestimate the 
total cluster mass, are unlikely to push the limit up to 
£Im = 1 ( Evrard 1997 and references cited therein ). 
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The region of parameter space which is allowed by 
these observational constraints for a flat model is then: 
0.5 < h < 0.7, 0.3 < n M < 0.6. These are also the limits 
that for MLM models the cosmological parameters must 
independently satisfy in order to achieve consistency with 
the set of linear clustering data previously analysed. 

We think that this is a notable feature of MLM mod- 
els and one of the most important results of this paper. 
We therL summarize our conclusions by saying that MLM 
models El appear to be a promising class of cosmological 
dark matter models. Our linear analysis shows that con- 
sistency for the cosmological parameters is achieved over 
a wide range of observational data. 
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